Feshbach resonances in an optical lattice 



in 
o 
o 

(N 
5h . 

< 
(N 



o 



CO 



I 

C 

O 

o 



(N 
> 

ON 

o 
o 

C3 



-a 
c 

o 
o 



X 
S3 



D.B.M. Dickerscheid/'E U. Al Khawaja, 2 D. van Oosten, 1 and H.T.C. Stoof 1 - 3 

institute for Theoretical Physics, Utrecht University, 
Leuvenlaan 4, 3584 CE Utrecht, The Netherlands 
8 Physics Department, United Arab Emirates University, P.O. Box 17551, Al-Ain, United Arab Emirates 
3 Department of Physics and Astronomy, SUNY, Stony Brook, NY 11794-3800, USA 

(Dated: February 2, 2008) 

We present the theory for ultracold atomic gases in an optical lattice near a Feshbach resonance. 
In the single-band approximation the theory describes atoms and molecules which can both tunnel 
through the lattice. Moreover, an avoided crossing between the two-atom and the molecular states 
occurs at every site. We determine the microscopic parameters of the generalized Hubbard model 
that describes this physics, using the experimentally known parameters of the Feshbach resonance 
in the absence of the optical lattice. As an application we also calculate the zero-temperature phase 
diagram of an atomic Bose gas in an optical lattice. 



I. INTRODUCTION. 

In the last few years there has been much excitement in 
the field of ultracold atomic gases. To a large extent this 
is due to two new experimental developments. The first is 
the use of so-called Feshbach resonances in the collision of 
two atoms, and the second is the use of an optical lattice. 
Both developments have led to an unprecedented cont- 
role over the physically relevant parameters of the atomic 
gas that can be used to explore new strongly-correlated 
regions of its phase diagram. In this paper we propose 
to combine these two developments and study an atomic 
gas in an optical lattice near a Feshbach resonance. 

A more specific motivation for studying Feshbach res- 
onances in an optical lattice is that recently it has been 
shown that in an atomic Bose gas near a Feshbach reso- 
nance a quantum phase transition occurs between a phase 
with only a molecular condensate (MC) and a phase with 
both an atomic and a molecular condensate (AC+MC) 
0, 0- The experimental observation of this quantum 
Ising transition is, however, complicated by the fact that 
in a harmonic trap the fast vibrational relaxation of Fes- 
hbach molecules consisting of two bosonic atoms appears 
to prevent the creation of a molecular condensate in that 
case |3- In an optical lattice with a low filling fraction 
molecule-molecule and atom-molecule collisions can es- 
sentially be neglected and we expect this problem to be 
much less severe. 

Having this particular application in mind, we from 
now on focus on atomic Bose gases. However, our re- 
sults can be immediately generalized to the case of a 
two-component Fermi gas or even an atomic Bose-Fermi 
mixture in an optical lattice. Moreover, we consider only 
such low filling fractions that it is justified to neglect the 
possibility of having three or more atoms per lattice site. 
The reason for this restriction is that in this case we have 
at most two atoms per site and the effect of the resonant 



interactions between the atoms can be incorporated into 
the theory exactly. The latter was shown previously to be 
very important for arriving at a quantitatively accurate 
description of a harmonically trapped atomic gas near a 
Feshbach resonance Q . How this can be achieved also in 
an optical lattice is discussed next. 



II. GENERALIZED HUBBARD MODEL. 

We consider the experimentally most interesting case 
of a deep optical lattice in which the on-site potential is, 
for low energies, well approximated by an isotropic har- 
monic potential with energy splitting hu> and the tun- 
neling energy t a for atoms between sites obeys i a <C fiio. 
For two atoms on a single site the two-channel Feshbach 
problem in the relative coordinate, after splitting off the 
center-of-mass motion, is then given by the Schrodinger 
equation 
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Here the noninteracting atomic Hamiltonian is Hq = 
— fi 2 V 2 /m + muj 2 r 2 /4. The bare detuning is denoted by 
(5b, r is the relative coordinate between the atoms and 
m is the atomic mass. The nonresonant or background 
atom-atom interaction is V aa and the atom-molecule cou- 
pling is denoted by V am . In first instance only the relative 
part is relevant, since only this part is affected by the in- 
teractions between the atoms. The center-of-mass part 
determines the tunneling. From Eq. Q we obtain the 
following equation for the molecules 
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where \ip m ) is the bare molecular wavefunction. Note 
that in the above we have implicitcly taken the extend 
of this wavefunction to be so small that its energy is 
not affected by the optical lattice, which is well justi- 
fied in practice. Because for most atoms we also have 
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that \V aa \ <C huj, we can neglect in the atomic propa- 
gator U aa compared to H . Moreover, the eigenstates 
\4> n ) of Hq with energy E n — (2n + 3/2)hu> that are rel- 
evant for an s-wave Feshbach resonance, can be writ- 
ten in terms of the generalized Laguerre polynomials 

as (l#„) = e-- 2 /4 i 2 i l/2 (r2/2z2)/(2 ^2 ) 3/4 [L l/2 (0)] l/2^ 

Here I — \J%/muj is the harmonic oscillator length. Us- 
ing these states Eq. (jSJ can be rewritten as 
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Using also the usual pseudopotential approximation, we 
have that (r| VamlV'm) = \/2gS(r), where the atom- 
molecule coupling g = fiyj2ira\ }S ABA^,/m depends on 
the background scattering length <2b g , the width of the 
resonance AB, and the difference in magnetic moments 
Aji of the relevant Feshbach resonance 4j. From this we 
then find that the energy of the molecules obeys 
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The function G(E) is the ratio of two gamma functions 
G{E) = T(-E/2hu + 3/4)/T(-E/2hu + 1/4). The di- 
vergence in Eq. Q, which was first obtained by Busch el 
al. in the context of a single-channel problem |(| , can be 
dealt with by using the following renormalisation proce- 
dure. The right-hand side of Eq. J3j can be interpreted 
as the selfenergy of the molecules KS m (E). The diver- 
gence in the selfenergy is energy-independent and is re- 
lated to an ultraviolet divergence that comes about be- 
cause we have used pseudopotentials. To deal with this 
divergence we have to use the renormalized detuning in- 
stead of the bare detuning. The former is defined as 
5 = 5b — linvio mg 2 /2ir?i 2 r, where 5 = A[x(B — Bo) is de- 
termined by the experimental value of the magnetic field 
Bq at resonance. Note that, as expected, the required 
subtraction is exactly equal to the one needed in the ab- 
sence of the optical lattice. In the latter case we have to 
subtract 2g 2 J dk m/7i 2 k 2 (27r) 3 HH, which can be in- 
terpreted as 6 — 5b — lim r jo 2g 2 / dk e lk ' r m/?i 2 k 2 (27r) 3 . 
In this manner we obtain the relative energy levels of 
the dressed molecules as a function of the experimental 
detuning that is shown in Fig. ^ 

From this figure we see that for very negative detuning 
the molecular state lies below the ground-state of the on- 
site microtrap and the bound-state energy is well approx- 
imated by the detuning. As it approaches the ground- 
state level of the trap there is an avoided crossing and 
as a result the lowest trap state is shifted upward. If the 
avoided crossings between the molecular level and sub- 
sequent trap states do not strongly overlap, the system 
can be well described by considering only the lowest trap 
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FIG. 1: The relative energy levels of the atom-molecule sys- 
tem as a function of the detuning S. This figure was calculated 
for g 2 /V2nl 3 (huj) 2 = 0.1 



state. The overlap between the avoided crossings is de- 
termined by the strength of the atom-molecule coupling 
and can be neglected if g 2 / ^/2tt1 3 (Huj) 2 <C 1. In this pa- 
per we restrict ourselves to a single-band approximation, 
although the generalization to the multi-band situation 
is straightforward. This means that we only take into 
account the wavefunctions of the molecular state and the 
ground state of the on-site microtrap. In that case only 
two energy levels are of importance when there are two 
atoms on a lattice site. We denote these levels by ej and 
e| and their behaviour as a function of detuning is shown 
in Fig. 

The effective atom-molecule coupling in the optical lat- 
tice is given by g> = g(J cfa# (x)| 4 )V2 = s/(2^ 2 ) 3/4 , 
where i/ , o( x ) is the Wannier function in the lowest band 
of the optical lattice. The effective atom-atom interac- 
tion is now given by U c s = U\> s — 2(g') 2 /(5 — 3fuv/2), 
where the background on-site interaction strength £/b g = 
(47rabgfi. 2 /m) J dx|?/>o(x)| 4 = \j2/ithu> (a^/l). It is in- 
teresting to note that in order for the single-band approx- 
imation to be valid we do not need to have that U e s -C hu 
because the on-site two-atom problem has been solved ex- 
actly. In Fig. we also show a close-up of the avoided 
crossing and the wavefunction renormalisation factors Z a 
that give the amplitude of the closed channel part of the 
molecules in the state \4> a ). Explicitely, we thus have in 
the single-band approximation that 
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The probability Z a is determined by the selfenergy 
of the molecules through the relation Z a = 1/(1 — 
dhY, m {E)/dE) 0. Note that in Fig. El the probabil- 
ity already shows the effect of the avoided crossing 
at a detuning of about Zhuj. As long as the single-band 
approximation is valid this will, however, not affect any 
of the results because the two-atom state that is involved 
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in this avoided crossing will not be populated. 

Combining the above we thus find a generalized Hub- 



bard Hamiltonian that is given by 
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Here t a and £ m are the tunneling amplitudes for the 
atoms and the molecules, respectively, and denotes 
a sum over nearest neighbours. The operators a|, a, 
correspond to the creation and annihilation operators of 
a single atom at site i respectively. The operators b\ a , 

b i correspond to the creation and annihilation opera- 
tors of the dressed molecules at site i respectively. Also 
e a = 3Tluj/2 is the on-site energy of a single atom. In the 
tight-binding limit the hopping amplitudes can be con- 
veniently expressed in terms of the lattice parameters as 
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Here A is the wavelength of the light used to create the op- 
tical lattice and lm_V2 = la = I- Note that, as expected, 
we have that t m oc t^/huj -C t a . Note also that our 
harmonic approximation to the on-site potential in prin- 
ciple slightly underestimates the hopping parameter. A 
more accurate determination of these parameters would 
involve the calculation of the appropriate Wannier func- 
tions. The chemical potential fj, is added because we per- 
form our next calculations in the grand-canonical ensem- 
ble. 



III. PHASE DIAGRAM. 

To find the mean-field phase diagram of a Bose gas 
in an optical lattice, we consider at sufficiently negative 
detuning the phase with only a Bose-Einstein conden- 
sate of molecules and perform a quadratic expansion of 
the Hamiltonian in the fluctuations of the molecular an- 
nihilation operator 6k, o- around the nonzero expectation 
value (&k, <r) = v /nmc ^ k -°^ cr .l- The effective Hamiltonian 
is then diagonalized by a Bogoliubov transformation and 
from the result we determine the equation of state of the 
gas as a function of the detuning 8 and the temperature 
T = l/kg[3. For the equation of state for the total filling 
fraction we find (cf. Ref. Q) n = n a + 2 with the 

molecular filling fractions obeying 




FIG. 2: Details of the physical content of our theory. We 
show the avoided crossing between the molecular level and 
the lowest two-atom trap state. The inset shows the prob- 
ability Z a as a function of the detuning 8. This figure was 
calculated for g 2 / y/2nl 3 (hoj) 2 = 0.1 Note that the center of 
mass contribution to the energy has been taken into account 
here. 



— ^mc 



±y I 

AT L^i e 0hu Kl _ i ' 



and the atomic filling fraction 
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Moreover, we have that N s is the total number of sites 
in the lattice, e£ = -2i a ff j=1 cos (k 3 -A/2) + e a , e£ CT = 
-2t m X^ =:l cos(fcjA/2) + £(Tj an d hu K<J = e£ ff + e m is 
the molecular dispersion. Likewise we find that hui^ = 
K e k — £ m/2) 2 - Ag^Z^n^c} 1 / 2 is the atomic Bogoliubov 
dispersion with e m = — zt m equal to twice the chemical 
potential and z is the number of nearest neighbours. 
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FIG. 3: Zero temperature phase diagram as a function of 
the filling fraction per site and the detuning 8 in units of hui. 
The different curves that separate the MC and the AC+MC 
phases correspond to values of g' /Ttw = 0.10 (full curve) and 
g' /hco = 0.12 (dashed curve) respectively. In both cases we 
have taken u> to be 10 4 rad/s. 

The critical temperature for the Bose-Einstein con- 
densation of the molecules follows from the condition 
n mc — 0. The location of the Ising quantum phase 
transitions follows from the zero-momentum instability 
in the atomic Bogoliubov dispersion when the detuning 
e m = — 4</y / Zj, l /rt mc + 2e a — 2zt a . In Fig. we show the 
results for this condition as a function of the total filling 
fraction and detuning. Note that in the limit of vanish- 
ing density the quantum critical point is determined by 
the ideal gas condition for Bose-Einstein condensation, 
i.e., [i = e m /2 — e a — zt a . From this condition it follows 
that for low enough filling fractions the location of the 
quantum phase transition shifts to higher detuning with 
increasing strength of the atom-molecule coupling. On 
the other hand at large negative detuning a larger value 
of the atom-molecule coupling implies a larger quantum 
depletion and hence a smaller molecular condensate frac- 



tion. This effect shifts the Ising transition to lower de- 
tuning. 

For completeness we would like to point out that at 
n = 1 the phase diagram can also contain a Mott- 
insulator phase |9(. This phase can occur at sufficiently 
large positive detuning such that U c s/zt a > 3 + 2\/2 
0, 13 • I n contrast to the quantum Ising transition, this 
transition has already been observed experimentally by 
Greiner et ol. flfjl after the theoretical prediction by 
Jaksch et ol. 11]. Its existence does not rely on the 
presence of the Feshbach resonance and we, therefore, 
have not included it in the phase diagram in Fig. |21 It 
is important to realize that this Mott insulator can only 
exist for repulsive interactions between the atoms, which 
requires (7b g to be positive. 



IV. CONCLUSIONS AND DISCUSSION. 

In summary, we have shown how to determine the mi- 
croscopic parameters of the generalized Hubbard model 
in Eq.@ that describes the physics of resonantly- 
interacting atoms in an optical lattice, using the experi- 
mentally known parameters of the Feshbach resonance in 
the absence of the optical lattice. As an application we 
also calculated the zero-temperature phase diagram of an 
atomic Bose gas in an optical lattice in the single-band 
approximation. By using an optical lattice one can sup- 
press three-body recombination processes that lead to a 
fast decay of the molecular condensate. 

For the single-band approximation to be valid the 
atom-molecule coupling constant g has to be small 
enough such that the avoided crossings between subse- 
quent bands do not overlap with each other. In some 
cases, however, this coupling constant can be too large 
for realistic conditions and the single-band approxima- 
tion will not hold anymore. In those cases we have to in- 
clude higher- lying two- atom states. To avoid this compli- 
cation the atom-molecule coupling g can be made smaller 
by using a more narrow Feshbach resonance or by us- 
ing two-photon Raman transitions to convert atoms into 
molecules 12]. Inclusion of higher- lying two-atom states 
is, however, easily achieved in our theory by adding more 
atomic and molecular states into the generalized Hub- 
bard model. In principle, we have to add several dressed 
molecular states for each additional atomic band 

that is required for a sufficiently accurate description of 
the atomic gas in the optical lattice. More precisely, for 
M atomic bands we need to include M 2 +l dressed molec- 
ular states into the theory. 
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